DIBARYON CONDENSATE IN NUCLEAR 



oo 



0^ 



X 



MATTER AND NEUTRON STARS: EXACT 
ANALYSIS IN ONE-DIMENSIONAL MODELS 

A. J. Buchmann and Amand Faessler 



^ ■ Institiit filr Theoretische Physik, Universitdt Tubingen 

Auf der Morgenstelle 14, D-72076 Tubingen, Germany 



M. I. Krivoruchenko 



; Institute for Theoretical and Experimental Physics 

O ; B.Cheremushkinskaya 25, 117259 Moscow, Russia 



Abstract 



Phase transitions of nuclear matter to the quark-gluon plasma with sub- 
O ' 

^ I sequent restoration of chiral symmetry have been widely discussed in the lit- 



erature. We investigate the possibility for occurrence of dense nuclear matter 
with a dibaryon Bose-Einstein condensate as an intermediate state below the 
quark-gluon phase transition. An exact analysis of this state of matter is pre- 
sented in a one-dimensional model. The analysis is based on a reduction of the 
quantization rules for the A^-body problem to coupled algebraic transcen- 
dental equations. We observe that when the Fermi momentum approaches the 
resonance momentum, the one-particle distribution function increases near the 
Fermi surface. When the Fermi momentum is increased beyond the resonance 
momentum, the equation of state becomes softer. The observed behavior can 
be interpreted in terms of formation of a Bose-Einstein condensate of two- 
fermion resonances (e.g. dibaryons). In cold nuclear matter, it should occur 
at 2{m]\i +ef) > w-d where m^r and rriD are the nucleon and dibaryon masses 
and ffi? is the nucleon Fermi energy. 



1 INTRODUCTION 



Exactly solvable models are of considerable interest. They provide important tests 
of different approximation schemes. These approximation methods may then be 



applied with increased confidence to more complex cases where analytic solutions are 
impossible. In fact, one-dimensional models very often provide the only possibility to 
gain analytic insights into the behavior of an interacting many-body system that go 
beyond perturbation theory. In this sense, exactly solvable low-dimensional models 
may provide valuable guidance in finding the proper dynamical description of more 
complicated systems. 

The recent discovery |T]] of Bose-Einstein condensation in a dilute vapor of 
rubidium-87 atoms created new interest in the phenomenon of Bose-Einstein con- 
densation and in exactly solvable models of interacting bosons. Already some time 
ago [0 it was noted that the energy spectrum of a one- dimensional Bose-gas of 
nonpenetrable particles is identical to that of an ideal Fermi gas. In other words, 
there is an exact one-to-one correspondence between bosons which interact via an 
"infinite wall" two-body potential, and a system of noninteracting fermions. Later 
P, Q an exact analytic solution for a one-dimensional Bose-gas interacting via a 
delta-function potential V{x) = a6{x), where a is an arbitrary positive parameter 
was obtained. It was shown that the model of ref. [2] is a special case of the model 
in refs. P, since in the limit a ^ oo particles do not penetrate each other. For 
a < oo the qualitative features of the energy spectrum remain unaltered. A class 
of exact solutions three particles with different masses interacting through a finite- 
strength delta-function potential was found ||^. Thermodynamic properties of an 
interacting Bose-gas with the potential V{x) = a6{x) were also discussed [^]. It was 
demonstrated that at finite temperatures the thermodynamic functions do not show 
a non-analytic behavior. In ref. the possibility of superfluidity at zero temper- 
ature is discussed. In ref. it was shown that long-range correlations inherent in 
the models of refs. ^ decrease as some power of 1/x at T = 0. One-dimensional 
models have also led to a deeper understanding of critical phenomena [^. Exact 
results for the quantization of the Toda lattice are reviewed in ref. [10]. 

There is an extensive literature devoted to dibaryon resonances. Dibaryons are 
predicted by most low-energy models of quantum chromodynamics (QCD) [11-14]. 
Some candidates are reported to be seen in experiments |]T3|, The appear- 

ance of Bose-Einstein condensation of two-fermion resonances in nuclear matter 
softens the equation of state of nuclear matter and decreases the upper limits for 
masses of neutron stars |]T^ . In the framework of the Walecka model, heterogenous 
nucleon-dibaryon matter is discussed in ref. |TH]. In the present paper, we focus on 
one-dimensional Fermi-systems with a Bose-resonance in the two-fermion channel. 
These one-dimensional models might have important implications for the behavior 
of dibaryon resonances in the nuclear medium. 

Previous results [2-8] were obtained for a 5-function interaction potential. Here, 



we show that one-dimensional models admit exact solutions for other types of zero- 
range singular potentials. We classify eigenfunctions of the A'"-body Hamiltonian, 
establish quantization rules, find eigenvalues of the A'"-body Hamiltonian without 
recourse to perturbation theory. In the thermodynamic limit one can exactly (nu- 
merically) calculate dispersion laws for elementary excitations, the ground state 
energy, and the equation of state. 

The outline of the paper is as follows. We start with discussing the simpler case 
of a system of bosons interacting through an arbitrary finite-range potential V{x) 
{V{x) 7^ for I a; |< a and V{x) = for | x |> a) generating a nontrivial two-body 
scattering S'-matrix. We then pass to the limit V{x) oo for | a; |< a and a ^ 0. 
In this limit, there is a wide class of nontrivial finite two-body S'-matrices. In sect. 
2, we establish the properties of the coefficients entering the plane wave expansion 
of the A^-boson wave function. In particular, we study their properties with respect 
to permutations of the particle quasi-momenta and establish relations between these 
coefficients. Then we discuss the symmetry of the wave function under permutations 
of the arguments, periodic boundary conditions, and matching conditions for the 
wave function. In Sect. 3, we derive the generalized Bohr-Sommerfeld quantization 
rules which reduce the quantization problem to the problem of finding solutions 
to N coupled algebraic transcendental equations. In the limit a — > 0, the energy 
eigenvalues are completely determined by the two-particle scattering phase shifts. In 
Sect. 4, we discuss the thermodynamic limit A?" — > oo, L — > oo, N/L — constant. A 
linear integral equation is derived for the distribution function of the quasi-momenta 
of the particles. The distribution function is completely determined by the scattering 
phase shifts of the particles. The elementary excitations in the Bose-system are 
classified and their dispersion laws are estabhshed. In Sect. 5, the results of Sects. 
2 and 3 are extended to special classes of exactly solvable Fermi-systems. In Sect. 
6, we discuss the properties of a system of fermions interacting via potentials which 
allow for a resonance in the two-body S'-matrix. In Sect. 7 the problem of stability 
of self- gravitating objects (neutron stars) made up of fermions with a resonance 
interaction is discussed. In one dimension all objects of such a kind are stable in 
general. In three-dimensional space narrow resonances cause instability of massive 
neutron stars. Finally, in sect. 8 we summarize and discuss the results. 



2 EXPANSION COEFFICIENTS OF THE WAVE 
FUNCTION 



We search for solutions of the iV-body Schrodinger equation in one dimension 

N a2 

(E(lf^ + E - -.^n) = E<f{xu xn) (2.1) 

in the interval [0, L]. Here, the particles with mass m are assumed to be bosons, so 
the wave function \E'(xi, xat) is symmetric under permutation of any pair of its 
arguments. The potential V{xi — xj) is assumed to vanish for \xi — Xj\ > a. Finally, 
we pass to the limit V{x) oo for | x |< a and a ^ 0. In this limit there exists a 
wide class of nontrivial two-body ^-matrices. 

The problem considered here is quite similar to the one- dimensional Heisenberg 
model of ferromagnetism with nearest neighbor interactions, first solved by Bethe 

Following refs. |@, ^, we impose periodic boundary conditions for the wave 
function 

^'(a;i, Xj = 0, ...,xn) = ^(a^i, xj = L, xn) (2.2) 

for any j. 

One can verify that the function 

N 

X{xi, Xn) = exp{i ^ kjXj) (2.3) 

j 

satisfies the Schrodinger equation, if xj + a < Xj+i, i.e. when V{xi — Xj) = for any 
pair of the arguments. The energy equals 

N h.2 

E = y-^. (2.4) 
^ 2m ^ ' 

3 

It is evident that any function which differs from Eq. (2.3) by a permutation of 
the particle quasi-momenta (fci, /ca?) — > (/cq^, fc^j^) satisfies the same equation 
and has the same energy. The problem reduces therefore to (i) the determination of 
the weights of all components {ka^, ka^f) of the wave function "^{xi, ...,xn), (ii) 
matching the wave functions in the different regions of integrability ("A-regions") 
where the exact wave function can be represented as a superposition of plane waves 
(2.3), (iii) symmetrization of the expression (2.3) with respect to the arguments, (iv) 
taking into account periodic boundary conditions, and finally (v) derivation of the 
multidimensional analog of the Bohr-Sommerfeld quantization rules for the particle 
quasi-momenta ki. 



Let us order the arguments of the wave function in increasing sequence, e.g. 
X5 < Xi < X3 < ... < X12. The numbers (5, 1,3, 12) constitute a set (71, ...,7jv). 
Therefore, 7^ is the number of the argument that occupies the i-th place in the 
above ordered sequence. Each A-region can be brought into correspondence with 
a set (7i,...,7jv) which is a permutation of the numbers (1,...,A^). There are N\ 
regions in which solutions can be represented in form of plane waves. Each 74-region 
is fixed by a set of inequalities 

x^. + a < x^.^^ (2.5) 

for J = 1, iV — 1, so that V{xi — Xj) = for any pair of the arguments. The solutions 
of the Schrodinger equation in the region (71, 7iv) can be written in the form 

N 

ai...ajv j 

The sum is taken over the A^! permutations of the (fci, /cjv)- It is assumed that 
the /c's are all different. It will be shown below that this is a general case. 

The wave function "i/^xi, ...,xn) contains for finite-range potentials terms with 
different sets of the particle quasi-momenta {ki}. In the limit a ^ 0, however, 
only one unique set {k^} survives. The requirement a — > is necessary to ensure 
completeness of the plane wave expansion in the ^-regions. 

Now, we consider the periodic boundary conditions. Suppose the set of argu- 
ments of the wave function on the left hand side of Eq. (2.2) belongs to the region 
(71, 7jv), then 71 = j due to the condition xj — 0. The set of arguments of 
the wave function on the right hand side of the Eq. (2.2) belongs to the region 
(Ai, Xn). 

It is clear that Xi — + 1 ior i < N, Xn — ^1 hy virtue of xj — L. The terms 
of the sum (2.6) are all hnearly independent. Therefore, the periodic boundary 
condition (2.2) can be unambiguously projected to the expansion coefficients 

C2^::Z:!^Z = C2tzi% eM^ka^L). (2.7) 

We illustrate this with an example. Let N = 2, j = 1 then 

*(0, X2) = CI2 exp{ik2X2) + C21 exp{ikiX2), 
^(L, X2) — C12 exp(ikiX2 + ik2L) + C|J exp{ik2X2 + ikiL). 

From the condition ^(0,^2) = ^(L, X2) follows C|| = C|| exp(iA;iL) and CH — 
Cll exp(i/c2l/). 

Under exchange of the coordinates x^^ ^ x^^ one region of integrability trans- 
forms to another one. The symmetry conditions for the wave functions under these 



permutations can be formulated in terms of the expansion coefficients in the follow- 
ing way 

Let us illustrate this again with an example. For N = 2, xi < X2 we get 

^{xi, X2) = exp{ikiXi + ik2X2) + CH exjp{ik2Xi + ik-iX2), 
^{x2i x\) — Cll exp(zA;iXi + ik2X2) + C|| exp(iA;2a^i + ikiX2). 

The condition '^{xi,X2) = '^{x2,Xi) implies C^l = 0^2, C^f = C^l- 

We consider the matching conditions for the wave functions in different A-regions. 
Consider first the two regions r± which can be obtained from each other by a 
permutation of the arguments x^^ and x^r such that s ±1 — r. In the sequence of 
increasing arguments of the wave function, the values x^^ and x^r occupy neighboring 
places. In the region r_|_, x^^ is on the right, i.e. s + 1 — r, while in the region r_, 
x^r is on the left, i.e. s — 1 — r. Let ^1 = x^s, — ^-yr- We now join the regions r+ 
and r_ and remove the restriction |^i — ^2] > o- In the region obtained (denoted 
herewith by F) , we are looking for solutions of the Schrodinger equation in the form 

N 

exp(i^%x,)x(6,6)- (2.9) 
3 

The sum is extended over j ^ s,r. Let gi = kas, ?2 = kar- In the region F, the 
function x(^i,^2) satisfies the equation 

(^ + §1 + ^(^1 - ^2))x(6, 6) = E'xiCi, 6) (2.10) 

2 2 

with E' = + 1^. The interval |^i — ^2! < ct, in which the potential is different 
from zero, is interesting only for supplying the correct matching conditions for the 
wave function for ^2 — > and ^2 — < In terms of the total and relative 
quasi-momenta of the particles, K — qi + q2, k — {q2 — ?i)/2, the center-of-mass 
coordinate X — {^1 + ^2)/2, and relative coordinate x — ^2 — ^i, the value x(^i, ^2) 
can be written for x < —a in the form 

exp(igi^i + ^^2^2) = ew{iKX + ikx). (2.11) 

If the incoming plane wave exp{ikx) is part of the solution, the outgoing plane wave 
exp{—ikx) exists as well: 

exp{iKX — ikx) — exp(igf2^i + ^^1^2). (2.12) 



These two waves differ from each other by permutation of the particle quasi-momenta 
gi and q2 only, and we do not get any additional solutions apart from the class of 



solutions of Eq. (2.6). Given that the expansion coefficients of the plane waves for 
X < —a (in the region r_) are known, one can reconstruct the expansion coefficients 
for X > a (in the region r+). This is a standard problem in scattering theory. 

There exist two linearly independent solutions which can be taken to be sym- 
metric and antisymmetric under the substitution x —x (^i <-> ^2): 

... X fr^yJ ^""'^ + S+{k)e-"'^; X < -a (region r_) 

X+(Ci,6) =exp(^i^X) <^ / • ^2.13) 

[ S+{k)e^''^ + e x > a (region r+) 

, I e^^^ — S-{k)e~'^^^: x < —a (region r_) 
X_(6,6)=exp(iKX) ^ / (2.14) 

[ S^_(A;)e*'^^ - e x > a (region T+j. 

The expansion coefficients in two neighboring regions of integrability are related by 
the scattering matrices S±{k). The scattering problem can be formulated on the 
half-axis x E (—00, 0] with the boundary conditions X+{^: 0)' — (symmetric case) 
and x_^(X, 0) = (antisymmetric case), or, equivalently, on the half-axis [0,-|-oo). 
The current density vanishes for x = for symmetric and antisymmetric wave 
functions. Prom particle number conservation it follows that the absolute values of 
the (S'-matrices are equal to unity on the real /c-axis. The following properties of the 
S -matrix hold true in the whole complex A;-plane: 

S±{k) = s±i-k)-\ 
s±{k) = {s±{k*)-^y. 

As usual, bound states are described by poles on the upper imaginary half-axis, 
virtual levels are described by poles on the lower imaginary half-axis. Poles in the 
lower half-plane of the complex /c-plane correspond to resonances. Por wave functions 
of the general form, x{^j ^) = C^x+i^^ ^) +C'_X-(^; x), the expansion coefficients 
for incoming e*^^ and outgoing e"*'^^ waves in the region r_ {x < —a) are connected 
unambiguously to the expansion coefficients in the region P_|_ {x > a). The relation 
is expressed through the scattering matrices S±{k) that can be obtained by solving 
equation (2.10) on the interval |^2~Ci| < Ci- We assume that this problem is already 
solved and that the matrices S±{k) are known. In Eq. (2.13) the signs of the S±{k) 
are fixed by the convention that for free particles S±{k) = 1 and X+(X, x) oc cos{kx), 
X-{X,x) oc sm{kx). Note that due to the boundary condition at a; = , the value 
X-{^) — exp(— iirX)x_(X, x) can be interpreted as the radial part of the scattering 
wave function in three dimensions, and the value S-{k), respectively, as the yS-matrix 
corresponding to zero angular momentum (the centrifugal potential for Z 7^ cannot 
be included in V{x) since V{x) is a short range potential). 

Identifying the components of the wave function x(X, x) — C+x+(X, x)-|-C_x_(X, x) 



with the relevant terms in the expansion of Eq.(2.6), we obtain 

C^::];:i;:X = C+S+ + C.S. (region r+, wave e^^-)^ 
C]^::1^X:X =C+-C. (region r+ . wave e-''^-) , 
<^5::-2;^-::7.", = C+ + C_ (region r_, wave e^'-), 

Q;:::ra;::X = ^+'^+ - C-S_ (region r_,wave e-'^-) 



(2.15) 



where S± = S±{{ka^ — A;„J/2). This system of equations is overdetermined. Yet, it 
admits a consistent solution for matching the wave functions. From these equations 
we derive 

9—9 2 

/^7i---7r7s---7]V _ .t_^__+r^7i---7s7r---7iv j /^7i---7s7r.--7]V (2 



Combining conditions (2.8) and (2.16), we obtain 



/^7l---7s7r---7JV _ C ((I. _ U W9\/^7i...7s7r.--7Af (0 A7\ 

If we make one more permutation of the indices a^, a^, we obtain, by virtue of 
S+{{ka^ — /cq,J/2) Sj^{{kas ~ kar)/'^) = 1; the initial expression which completes this 
consistency check. 

3 GENERALIZED BOHR-SOMMERFELD 
QUANTIZATION RULES 

It follows from condition (2.17) that the expansion coefficients are all expressed 
unambiguously through Cl'"^. Using Eq.(2.8), the upper indices can be ordered In 
sequence (1,...,A^). The expansion coefficients, which can be obtained from each 
other by permutation of two neighboring indices, are related by Eq.(2.17). The 
arbitrary set of indices (cti, «Ar) can be obtained from the sequence (1, N) by 
permutations of neighboring indices. Thus, one can express C^;^'^^^ through Cl'"^. 
It remains to verify that different sequences of transpositions transforming (1, A^) 
to («!,..., ttjv) give the same result. Let Pi and P2 be two such permutations of 
the initial sequence (1, .....N) leading to the same final sequence (ai, Oat). It is 
evident that the permutation Pi x P2~^ describes an identical (trivial) transformation. 
In such a sequence, each pair of indices changes places an even number of times. 
Since S+{{ka^ — ka^)/2) S^{{ka^ — ka^)/2) — 1, the result of the transformation 
Pi X P2-' is the identity Cl;;;^ — C};;;^. Denoting the result of transformation Pi 

by 

cl;^l. = act:^ (3.1) 



and the result of transformation P2 by 

Clz^N = f3C'^;l. (3.2) 

we obtain a/3 = 1. Therefore, from both permutations Pi and P2 we obtain one 
and the same relation (3.1). The phase shift 6+{k) is defined by equation S^{k) = 
exp{2i6+{k)). With the help of eqs. (2.7), (2.8) and (2.17) we obtain the generalized 
Bohr-Sommerfeld quantization rule 

k,L + Y,26+{^ -) = 27rn,-. (3.3) 

1=1 2 

The sum is running over / 7^ j, nj are integer numbers. This equation can 
be interpreted in the following way. Going around the circle [0,L], the particle j 
scatters on each particle I ^ j, acquiring an additional phase 2^+( *'^~*^' ). When the 
particle j is back to the initial place, its phase turns out to be equal to the left hand 
side of Eq.( p.3| ), the additional term kjL is a result of the translation. The total 
phase must be an integer multiple of 27r, since the wave function is single valued. 

Let us now consider the case when one virtual level exists in the complex k 
-plane. In this case the S'-matrix 

S+{k) = {k-iko)/{k + iko), (3.4) 

where ko = ma > 0, corresponds to the delta-function potential V{x) = a6{x) 
(problem 2.47 in ref. 0). Respectively, 5+ (-00) = 0, 5+ (+00) = 2n, 6+{0) = vr, 
so that 5*4.(0) = —1. It is seen from eqs. (2. 13) that the symmetric wave function 
vanishes for /c = and ^^(O) = —1 and the particle quasi- momenta in the set 
{ki, k2, kjy) must be all distinct. Therefore, one obtains a constraint which re- 
sembles the Pauli principle even though we have started with a system of bosons. In 
the lowest energy state, the particle quasi-momenta (fci, k2, /cat) occupy the Fermi 
sphere. The equality ^^(O) = — 1, along with the exclusion principle, apparently, is 
valid for any odd number of virtual levels. 

For an even (zero) number of virtual levels one has ^^(O) = 1, and some quasi- 
momenta in the set {ki, k2, kjy) can coincide. However, for qi = q2, \x\ > a, Eq. 
(2.10) has the general solution x+(X, x) = (a + bx) exp{iKX) with 6 7^ 0. For this 
reason it is impossible to satisfy the periodic boundary condition. The states with 
b ^ correspond to zero-energy scattering states. The coefficient b vanishes for a 
special class of potentials having a representation of the form 

Vix) = -x"{x)/xix) (3.5) 
m 

with the wave function satisfying the condition x'(0) = x'{(^) = 0. From V{a) = 
it also follows that x"{(^) = 0. Solutions of the Schrodinger equation with a limited 



asymptotic behavior at infinity (i.e. 6 = 0) exist when a discrete level in the potential 
appears (problem 2.18 in ref. [0). Therefore, it is clear that the case 6 = is an 
exceptional one. Solutions of such a kind occur when we pass to a system in which 
bound states are formed. 

We thus consider repulsive potentials without discrete states in the energy spec- 
trum. For such potentials all quasi-momenta in the set {ki, k2, k^) must be 
distinct. 

4 THERMODYNAMIC LIMIT 

In the thermodynamic limit N oo , L ^ oo , and p = N/L a fixed value, the 
particles in the ground state occupy continuously the Fermi sphere. The sum over 
the particle quasi-momenta can be approximated by an integral over k with a weight 
function f{k) 

where Lf{k)dk/{2TT) is the number of states in the interval dk. The distribution 
function f{k) can be found from equation 

m = i+ r ^^m5'/--^) (4.1) 

where (5^(^^) is the derivative of the scattering phase shift with respect to the 
argument. This equation is obtained by rewriting Eq.(3.3) in the thermodynamic 
limit for the phase difference of the particles j + 1 and j taking into account that 
rij+i — rij = 1. Eq.(4.1) is a generalization of Eq.(3.12) in ref. [3] to arbitrary 
singular potentials. 

We consider the problem of finding the spectrum of elementary excitations. Let 
us consider two sets of particle quasi-momenta {ki, k2, /ca?) and {k[, k'2, k'^). In 
the first set the /c's occupy continuously the Fermi sphere. The second set of fc's is 
obtained from the first one by removing a particle from the Fermi surface with the 
quasi-momentum k^ = Pf and giving it an arbitrary quasi-momentum q = k'^ > pp 
OT q = k'^ < —pF- Because of the interaction, the particle quasi-momenta inside the 
Fermi sphere receive a shift kj — kj = A{kj)/L. Taking the difference between the 
two relations 

k,L + E 254^) + 25+(^S^) = 2™,, 
1=1 ^ ^ 

1=1 ^ ^ 



we obtain with the help of Eq.(4.1), 

A{k)f{k) = -254^) + 25+(^) + ^A{k')f{k')6'4'^). (4.2) 
The total momentum of the system equals 

/Pf dh' 
—A{k')f{k'). (4.3) 

The first term is the quasi-momentum of the excited particle. After removing the 
particle with quasi-momentum pp the momentum of the Fermi sphere changes by 
an amount —pp- The last term represents the change of the total momentum of the 
particles inside the Fermi sphere. The energy of the system equals 

e=f-^+r^-^ik')m. (4.4) 

With q ^ Pf the momentum P and the energy e of the excited state vanish. In order 
to find the dispersion law for elementary excitations, e(P), it is necessary to express 
q in terms of P with the help of Eq.(4.3) and substitute the resulting expression into 
Eq.(4.4). Note that the integral equations (4.1) and (4.2) have the same kernel 

R{k, k') = {2n)6{k - k') - 5;(^^) (4.5) 

and that R{k, k') is symmetric. Solutions of these equations can be represented in 
the form 

/Pf fjk' 
—R-\k,k'), (4.6) 
-pF 217 

mm = ^R'\k, k'){-26^{^) + 264^^)). (4.7) 

Integrating Eq. ( |4.7| ) over k and using the symmetry under permutation of the argu- 
ments of R~^{k, k'), we can rewrite Eq.(4.3) in the form 

fPF dk' j,^^!^^ ^k'-q^ ^ ^k'-pp, 

-pF 

In lowest order of the difference q — Pf 



P=,_p,+ r f^f^k')^-2SA'^)^2SA'^)y (4.8) 

J-PF 2tt 2 2 



P = {q-PF)f{PF), (4.9) 



/Pf rlh' h' 
^-R-\k\pF). (4.10) 
-PF m. 

In a similar way one can treat excitations which transfer a particle to the opposite 
side of the Fermi surface, that is for ki = —pp and q = k[ > pp or q = k[ < —pp. 
These excitations are described by the same formulae, if the substitution pp —pp 
is made. 



Other kinds of elementary excitations of "hole" type (non-Bogohubov excita- 
tions) are obtained by taking away a particle with quasi-momentum q inside of the 
Fermi sphere and placing it on the Fermi surface. Let us find the spectrum of such 
excitations. Consider two sets of the particle quasi-momenta (/ci, A;2, /cat) and 
{k[, k'^). In the first set the quasi-momenta occupy continuously the Fermi 
sphere, while the second set is obtained from the first one by removing a particle 
with a quasi-momentum k,n = q from the Fermi sphere and placing it on the Fermi 
surface at k!^ = Pf- By comparison of the relations 

k,L + E 25+(^) + 254^) = 2n„ 
1=1 ^ ^ 

+ E 2M^) + 2M^) = 2n, 

with the corresponding relations for "particle" type excitations, discussed before, 
we see that the corresponding equations can be obtained from the ones already 
discussed by the replacement Pf ^ q- In the above formulae the terms I — m are 
excluded from the summation. 

In the symmetric case, when the particle with quasi-momentum km — q^^ placed 
on the Fermi surface at k'^ — —pr, it is necessary to make the replacement Pf ^ 
-Pf- 

5 EXACTLY SOLVABLE MODELS FOR FERMI- 
SYSTEMS 

In full analogy to the Bose case, the problem of constructing the wave function of 
the Fermi-system can be solved when the spins of the fermions are all lined up in 
one direction. In such a case, the spin part of the wave function is symmetric, and 
the coordinate part is totally antisymmetric. Let us find how eqs.(2.7), (2.8), (2.17) 
and (3.3) should be modified. 

The periodic boundary condition (2.7), apparently, remains unaltered. 

Under permutations of the coordinates, the wave function is antisymmetric. The 
symmetry conditions yield the expansion coefficients satisfying the condition: 

/^7l...7s...7r-.7Ar _ _flll--yr-ls-lN 

In comparison to the Bose case, the right hand side acquires a minus sign. 

The matching conditions of the wave function at the boundaries of the differ- 
ent regions of integr ability are derived as for bosons. Relations (2.16) are valid. 



(5.1) 



Combining eqs.(2.16) and (5.1), we obtain 



^ai...asar-aN ~ '-^-yy'^'^r '^ocs) I ^)^ai...aras-aN- y^"^) 

The expansion coefficients are related to each other through the matrix S^{k) 
that describes scattering on the positive half-axis with the boundary condition 

x(o) = 0. 

The generalized Bohr-Sommerfeld quantization rules have a form analogous to 
Eq.(3.3) 

^ k - -k, 

kjL + ^2S4^-^)^2nn,. (5.3) 

1=1 ^ 

where the phase shift is defined by S^{k) = exp(2i5_(fc)). The summation is per- 
formed over I ^ j, rij are integer numbers. 

Therefore, there is a close analogy between the the behavior of Bose- and Fermi- 
systems with parallel spins. This analogy also exists in the thermodynamic limit. 
The results of Sect.4 are valid for Fermi-systems after the replacement S+{k) ^ 
S.{k). 

In the limit /c — > 0, the scattering phase shifts S+{k) and S-{k) have a differ- 
ent behavior. For smooth potentials, the continuity conditions for the logarithmic 
derivatives at x = a in the symmetric and antisymmetric cases have the form 

ktan{ka + 5+{k)) = K,+, (5.4) 

kcot{ka + 6-{k)) = K-, (5.5) 

A;+ and /c_ do not depend on the momentum k if /c << h_. We conclude that 
at small k the scattering phase shifts have the form 

5+(/c) = —ka + arctan(fi;+//c), (5.6) 

5-{k) — —ka -\- arccot(«;_/A;). (5.7) 
In the limit a — > and for arbitrary small but finite values of and , we obtain 

5+{k) 7re{k), (5.8) 

54k) 0. (5.9) 

In the weak coupling regime (the value of the k+ is small) the quantization conditions 
(2.17) for the symmetric case remain non-trivial, since the derivative of d^{k) is 
proportional to a delta-function. The integral kernel R{k, k') defined by Eq.(4.5) is 
determined by the difference of the delta-function and the derivative of S+{k). For 



Bose-systems with the potential V{x) — aS{x) , it is impossible to reproduce the 
result of perturbation theory analytically beyond first order of a [4] , because of the 
comphcated character of the kernel R{k,k'). In general it is, however, not difficult 
to perform a comparison with perturbation theory numerically. The non-trivial 
character of the weak coupling regime can be related to the non-analyticity in the 
coupling constant (in A if we write the potential as \V{x)), that should be present 
because of the instabihty of the system of bosons when the sign of the potential is 
changed. 

According to Eq.(5.9) the phase 5^{k) and its derivative should be set equal to 
zero. After that, the quantization conditions (5.3) take a very simple form kjL = 
21111 j. As a result, we deal with a Fermi-gas of non-interacting particles. The weak 
coupling regime is therefore trivial for fermions. 



6 BOSE-EINSTEIN CONDENSATION OF TWO- 
FERMION RESONANCES 



In this sect, we consider a problem of considerable physical interest. As in sect. 
5 consider a one-dimensional system of fermions interacting via a finite ranged po- 
tential V{x). Suppose there exists in the fermion-fermion channel a resonance with 
momentum ko = ki — ik2, where ki and k2 are real, positive numbers. One expects 
that after increasing the Fermi momentum kp above the value of the resonance ki, 
the creation of such resonances will be energetically favored by the system. The 
resonances can be treated as composite Bose particles. If their interaction is small, 
they are accumulated in the ground state like real bosons with zero total momen- 
tum. The Bose-Einstein condensation reveals itself by an increase of the distribution 
function f{k) of the particle quasi-momenta near the Fermi surface, since the res- 
onances are at rest and the fermion momenta are concentrated in the vicinity of 
k ^ ki ^ kp- The effect is well pronounced for a small width and disappears with 
increasing width of the resonance. 

Near the resonance, the (S-matrix can be parametrized in the Breit-Wigner form 



An S'-matrix of such a type corresponds not only to finite range potentials, but 
also to zero-range singular potentials. 
Let us consider a potential of the form 




(6.1) 



V{x) 



Vo9{a — x) + a6{x — a), 



(6.2) 



where Vq represents a positive value. We wish to find a solution of the scattering 
problem on the half-axis [0, +00). The wave function takes the form 

X(x)^{ ^ ^ ' (6.3) 

[ Csm{kx + d_{k)); a < x. 

Here, K — ^/k^ + Vq (in "atomic units" m — h — 1). The scattering phase is 
determined by the condition 

kcot{ka + S-{k)) = K cot{Ka) + a. (6.4) 

The ^'-matrix takes the form 

S4k) = exp(-2iH ^^"^;^"|+" + t (6-5) 
^ ^ ^ K cot{Ka) + a - ik ^ ' 

The (S'-matrix poles can be found from the equation 

tan(ira) = ., ^ . (6.6) 
ik — a 

In the limit a — > 00 , the (S-matrix poles cross the real axis at Ka — {mod tt 
). The lowest level is located at Ka — tt . We may keep the lowest resonance 
momentum ki — y^TT^/a^ — Vq fixed and pass to the limit Vq — > 00 and a — > 0. The 
zero-width resonance occurs then from a solution to the Schrodinger equation with 
the zero-range singular potential (6.2). The potential is determined by the limiting 
procedure: Vq 00, a ^ 0, /ci = ^tt^/o^ — Vq = constant. It generates on the real 
A;-axis exactly two poles (one zero-width resonance) in the 5* -matrix. The other 
poles are moved to infinity. 

Let now a be finite, but large, and Vq receives a correction AVq << Vq. The 
resonance acquires, first of all, a finite width. We require the AVq be such that the 
real part of the resonance momentum be equal to ki. The ^'-matrix pole is located 
at ko = ki — ik2. 

Now, we parametrize a = l/{^x^) and a = x^, so that 

Voy ^-,-kl (6.7) 
x^ 



In the hmit x — > 0, Eq.(6.6) gives 
AVo/n' = + ^ _ ^(6 - - 3x'kl) + ^(10 - IStt^ - 27x'kl) + ... , (6.8) 

X X oX o 

k^/TT^ = - 3xe' + (6 - tt" - x^kl)x^^' + ... . (6.9) 

As a result, we obtain a zero-range singular potential for which the 5* -matrix in the 
complex /c-plane has two poles only, which can be identified with a resonance. The 
resonance width is proportional to For narrow resonances ^ << 1. 



In Fig.l, we show results for the distribution function f{k) for different values 
of the Fermi momentum kp from 0.64 to 1.6 with a step 0.12 in a system where the 
two-fermion interaction is described by a S-matrix of the form (6.1) with ki = 1, 
/c2 = 0.05 (small width). The distribution function f{k) increases near the Fermi 
surface where the production of the resonances becomes energetically favorable. This 
effect can be interpreted in terms of a Bose-Einstein condensation of the resonances 
whose wave function are concentrated aX k ^ ki ^ kp- The value k2 measures 
the spread of the fermion momenta in the resonances. When the Fermi momentum 
is increased beyond the resonance momentum ki the distribution function shows a 
plateau centered at /c ~ fci. If all resonances are in the Bose-Einstein condensate, 
the size of the plateau would be of order k2- However, its size increases with the 
Fermi momentum. This can be interpreted as follows: It is known that in systems 
of interacting bosons there exists a fraction of bosons with nonvanishing velocities 
which are not in the condensate and which have a nontrivial momentum distribution 
pT |. These bosons contribute to the pressure. Their fraction increases with the total 
density of the bosons. 

The appearance of resonances in Fermi systems yields a softer equation of state. 
This means that the pressure p increases slower with the density n . The effect is 
shown in Fig. 2. The sudden rise oi p = p{n) appears at kp = ki where the particle 
density equals n = 0.36. 

For a broad resonance, the derivative of the scattering phase S' {^^) entering 
Eq.(4.2) is small, the distribution function f{k) is close to unity, and therfore the 
effect of the resonance on the equation of state is small, too. 



7 SELF-GRAVITATING OBJECTS 

Bose-Einstein condensation of narrow two-fermion resonances may drastically change 
the properties of fermion matter, producing physically interesting phenomena, for 
example, the instability of neutron stars. Before, studying this effect in more detail 
let us qualitatively describe the expected scenario. 

Narrow two-fermion resonances can be treated as Bose particles. If the central 
density of a neutron star exceeds a critical value, creation of these bosons with 
subsequent formation of a Bose-Einstein condensate becomes energetically favorable. 
The critical density for a Bose-Einstein condensate formation is determined by the 
mass of the resonances. In the ideal Bose-gas approximation, the chemical potential 
of the fermions is frozen at /i = m£)/2, where mo is the dibaryon mass. The 
dibaryons are at rest. Therefore, they do not collide with the boundary and do not 
contribute to the pressure. The pressure is determined by the fermions only. The 



number of fermions remains constant with increasing density, since the radius of the 
Fermi sphere is frozen, whereas the number of the resonances increases hnearly. The 
incompressibihty of matter vanishes. In this way, we give a quahtative explanation 
for the observed growth of the distribution function f{k) near the Fermi surface in 
Fig. 1 as well as the rapid change of p = p{n) in Fig. 2. 

Suppose that the short-range potential between fermions is such that a narrow 
two-fermion (dibaryon) resonance exists and that a Bose-Einstein condensate of 
these resonances is formed in the interior of a neutron star (in three dimensions) for 
r < ri. In the inner region, due to the formation of the resonances, the pressure 
remains constant, i.e. Vp = 0. Gauss's law implies 



J dS ■ V0(r) = 47rGM(r) 



where M(r) is mass of the substance inside of a sphere of the radius r and G the 
gravitational constant. We conclude that V0(r) 7^ 0. The Euler equation 

dv 

p-^ + p{v-V)v = -Vp-pV(j){r) (7.1) 

shows that in the static case (v = 0) the gradient of the pressure p is balanced by 
the gradient of the gravitational potential (f){r). However, if Vp = 0, static solutions 
are impossible and neutron stars are gravitationally unstable. 

Let us study this effect in a one-dimensional model of the neutron star, where all 
quantities can be calculated exactly, in order to test if some instability occurs due to 
the formation of the two-fermion resonances. We restrict ourselves to nonrelativistic 
Newtonian gravity. The static stable configurations of neutron stars are described 
by Euler 's equation with a vanishing left hand side 

-^-p(x)^ = (7.2) 

where p{x) is the pressure, p{x) = mn{x) the mass density, n{x) the number density, 
and (t){x) the gravitational potential. The center of the neutron star is placed at the 
origin of the coordinates. Gauss's law gives 

= 47rG'p(x). (7.3) 



dx'^ 

Note that in contrast to classical electrodynamics the right hand side of Eq.(|T] 
has a positive sign. 

We integrate Eq.(7.3) from —x to +x. Using the sjnnmetry of the potential (t){x) 
under the reflection x ^ —x {p{x) and p{x) are also symmetric functions) and the 
antisymmetry of the first derivative of (f){x), one gets 

= 27rGM(x), (7.4) 

dx 







where 

M(x) = j dxp{x) = 2 I dxp{x). (7.5) 
Substituting Eq.(7.4) into Eq.(7.2), we obtain 

dp 



dM 

and finally 



+ ttGM = 0, (7.6) 



M(x) = ^-^(p(0)-p(x)). (7.7) 

At the surface p{xs) = 0, and the total mass of the neutron star Mg = M{xs) is 
expressed unambiguously through the central pressure p{0). The neutron star radius 
Xs can be obtained from Eq.(7.7). It is sufficient to take the derivative of Eq.(7.7) 
with respect to x, divide both sides by 2p(x), and integrate the result over x. In 
this way one gets 

= -J— r^'^ (7.8) 

V2^Jo p{p)^p{0)-p 

Assuming that the equation of state n = n{p) is known, eqs.(7.7) and (7.8) 
allow to determine the neutron star radius function of the neutron star 

mass Msi or, equivalently, the total mass Ms as a function of the central pressure 



]9(0). The criterion for the gravitational stability of stars has the form (see ref. |22 
Eq.(10.1.4p)) 

where x{M) is a coordinate, such that inside of the interval [— x, x] the matter of 
mass M is contained: 

, , dM' 
.(M)=l ^. (7,10) 

Taking the derivative with respect to Ms and using Eq. (7.7) to express the pressure 
in terms of the mass, we obtain 

dx{M) ttGMs f^^ dM' dp 



dMs 2 Jo p2 

In the non-degenerate case dp/ dp = > (a^ is the velocity of sound), the deriva- 
tive dx{M) / dMs is negative definite and one-dimensional neutron stars are in gen- 
eral gravitationally stable. 

In three-dimensional space this is not necessarily the case. In order to investigate 
the reason for this difference, let us express the sum of the internal energy and the 
gravitational binding energy 

E = j edV + sign(2 - D)G j r^-^MdM (7.12) 



in terms of the average density p ~ M/r^ where D is the dimension of the space. 
For an equation of state of the form e = Ap"' we get 

E = AfP-^M + sign(2 - D)BM(2+^)/^p(^-2)/D (^^^g) 

where i? is a fixed positive constant. 

In the one-dimensional case D = 1, a minimum energy as a function of p is 
obtained for 7 > 1. In the ideal gas approximation for the Bose-Einstein condensate 

7 = 1, p ~ mun, where n is the number density, and so there is no stability. The 
maximum pressure determines the maximum mass of the neutron star 

M^ax = ^^Pmax(O). (7.14) 

However, if one goes beyond the ideal gas approximation, after the creation of the 
resonances the pressure still increases slowly with the density due to the interac- 
tions between the fermions and bosons as shown in Fig. 2. In the one-dimensional 
problem, even an arbitrarily slow growth of the pressure (7 = 1 -1-0.0. ..01) stabilizes 
the neutron star. 

In the three-dimensional case, the minimum of E exists for 7 > 4/3, and for the 
stabilization of the neutron star, an arbitrarily slow growth of pressure is insufficient. 
In order to make the neutron star stable, it is necessary to increase the average value 
of 7 above 4/3 (see [^). Consequently, in the presence of a dibaryon condensate 
a considerably stiffer equation of state is required to prevent the neutron star from 
collapsing. 

8 CONCLUSION AND DISCUSSIONS 

We have shown that the properties of one-dimensional Bose- and some Fermi-systems 
can be determined exactly without recourse to perturbation theory for a wide class 
of singular potentials. The solutions are expressible through two-particle scattering 
phase shifts. The Fermi character of the energy spectrum in one-dimensional Bose- 
systems is not specific to potentials of delta-function type. In any singular potential 
there are hole-type non-Bogoliubov branches of elementary excitations. We have 
derived the dispersion laws for these excitations. 

The one-dimensional Fermi system can be analyzed exactly if the fermion spins 
are lined up in one direction. The quantization rules have a similar form for Bose- 
and Fermi-systems. 

An exact analysis of Fermi-systems with a resonance in the two-fermion channel 
is given. In the model considered, the Pauli principle and the composite nature of 



the resonances are taken into account at the outset. We have observed an increase of 
the distribution function f{k) of fermions over the quasi- momenta k in the vicinity 
of the Fermi surface when the density is close to the critical density of resonance 
formation. We could interpret this behavior in terms of a Bose-Einstein condensation 
of two-fermion resonances. The formation of the resonances is accompanied by a 
softening of the equation of state. In the real (three-dimensional) world softening of 
the equation of state of nuclear matter caused by dibaryon resonances can produce 
instability of neutron stars. In a one-dimcnsional model of self-gravitating objects 
this effect does not exist. 

The existence of a dibaryon resonance will lead at higher densities to the occur- 
rence of a new state of nuclear matter. The dibaryons behave like bosons and can 
form a Bose-Einstein condensate in nuclear matter. With increasing baryon density, 
the pressure should only slightly increase with the density. In a perturbative pic- 
ture (which was not used here), the dibaryon condensation should occur at a Fermi 
energy ep determined by the condition 

2{mN + sp) > rriD 

where mjv and niD are the nucleon and dibaryon masses. 

It will be interesting to search for such a dibaryon condensate in heavy-ion col- 
lisions. In the center-of-mass frame of the condensate a large fraction of dibaryons 
has zero velocities. When the density decreases, dibaryons in the condensate decay 
into nucleons. If a channel D NN exists, experimentalists will observe monochro- 
matic nucleons with the energy mD/2 . Since the rest frame of the condensate is a 
priori unknown, it is necessary to look for a boost transformation along the beam 
momentum into a coordinate system in which a statistically significant excess of 
monochromatic nucleons exists. An excess of such events can be considered as a 
possible signature for the formation of a dibaryon condensate. 

It is also interesting to check astrophysical data for the presence of a dibaryon 
condensate in the interiors of massive neutron stars as well as for signatures of 
instability of neutron stars caused by dibaryons. 
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FIGURE CAPTIONS 
Fig.l The distribution function f{k) versus the quasi-momentum k < kp for 8 
different Fermi momenta kp ranging from 0.64 to 1.6 with a step 0.12 (in units 
m — h = 1). The two-body ^'-matrix has a pole corresponding to a resonance at 
k — ki — ik^ with ki — 1 and k2 — 0.05 (i.e. small width). The distribution function 
f{k) increases near the Fermi surface when kp approaches ki. This effect can be 
interpreted in terms of a Bose-Einstein condensation of the resonances, since the 
wave function of the resonance is concentrated at /c ~ /ci. 



Fig. 2 Pressure p versus fcrmion number density n. The change of the slope occurs 
ai n — 0.36 when kp = ki (in units m = h = 1). This is the point at which the 
two-fermion condensation starts. It results in a softening of the equation of state. 
In the ideal gas approximation, the pressure does not increase with the density 
after the condensation has set in. If interactions are included, the pressure is slowly 
increasing even after the condensation. 



